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ABSTRACT 

The High Resolution Fly's Eye (HiRes) experiment is an air fluorescence detector which, operating 
in stereo mode, has a typical angular resolution of 0.6° and is sensitive to cosmic rays with energies 
above 10 18 eV. HiRes is thus an excellent instrument for the study of the arrival directions of ultrahigh 
energy cosmic rays. We present the results of a search for anisotropies in the distribution of arrival 
directions on small scales (< 5°) and at the highest energies (> 10 19 eV). The search is based on 
data recorded between 1999 December and 2004 January, with a total of 271 events above 10 19 eV. 
No small-scale anisotropy is found, and the strongest clustering found in the HiRes stereo data is 
consistent at the 52 % level with the null hypothesis of isotropically distributed arrival directions. 
Subject headings: cosmic rays — acceleration of particles — large-scale structure of universe 



1. INTRODUCTION 

Identifying the sources of ultrahigh energy cosmic rays 
remains one of the central challenges in astrophysics. Af- 
ter three decades of systematic searches for the origin of 
these particles, source identification still remains elusive. 
Sky maps of cosmic ray arrival directions at all energies 
are generally isotropic, with no obvious source or source 
region standing out. 

A direct way to search for sources of ultrahigh energy 
cosmic rays is to analyze the distribution of their ar- 
rival directions for small-scale clustering. Any signifi- 
cant clustering in arrival directions could be evidence of 
nearby, compact sources, whereas the lack of clustering 
is consistent with models in which ultrahigh energy cos- 
mic ray sources are distributed at large distances from 
our Galaxy. 

Arrival directions do not necessarily point back to 
sources, as charged cosmic ray primaries suffer deflec- 
tions traveling through Galactic and intergalactic mag- 
netic fields. The strength and orientation of these fields 
is not well established, so the size and direction of the 
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deflection is difficult to ascertain. However, since the 
Larmor radius increases with energy, the possibility of 
observing small-scale anisotropy associated with cosmic 
rays pointing back to their origins is expected to grow. 

Indeed, small-scale clustering of cosmic ray arrival 
directions at the highest energies has been previously 
claimed. The AGASA (Akeno Giant Air Shower Array) 
experiment reported possible clusterin g in their sample 
of eve nts with energies above 4 TO 19 eV l)Havashida et al.l 
1996). The analysis has been updated severa l times 
l|Takeda et aDll999t l200lt iTeshima et al.l 1200^ . most 
recently reporting six clusters (five doublets and one 
triplet) in a sample of 59 events, where a cluster is de- 
fined as a set of events with angular separation less than 
2.5°. The chance pro bability of this signal was reported 
to be less than 10" 4 l|Teshima et al. 11200.^ . 

Given the potential importance of this result for our 
understanding of the origin of cosmic rays, it is crucial to 
test the claim that clustering is a feature of cosmic ray 
arrival directions with independent experimental data. 
Since 1999, the High Resolution Fly's Eye (HiRes) air flu- 
orescence experiment has been operating in stereo mode, 
collecting data of unprecedented quality on the arrival 
direction, energy, and composition of ultrahigh energy 
cosmic rays. In this Letter, we report results of a search 
for small-scale anisotropy in the arrival directions of ul- 
trahigh energy cosmic rays observed by the HiRes stereo 
detector between 1999 December and 2004 January. 

2. THE HIRES DETECTOR 

HiRes is an air fluorescence experiment with two sites 
(HiRes 1 & 2) at the US Army Dugway Proving Ground 
in the Utah desert (112° W longitude, 40° N latitude, ver- 
tical atmospheric depth 860g/cm 2 ). The two sites are 
separated by a distance of 12.6 km. 

Each of the two HiRes "eyes" comprises several tele- 
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scope units monitoring different parts of the night sky. 
With 22 (42) telescopes with 256 photomultiplier tubes 
each at the first (second) site, the full detector covers 
about 360° (336°) in azimuth and 3° - 16.5° (3° - 30°) 
in elevation above horizon. Each telescope consists of a 
mirror with an area of about 5 m 2 area for light collection 
and a cluster of photomultiplier tubes in the focal plane. 

A cosmic ray primary interacting in the upper atmo- 
sphere induces an extensive air shower which the detec- 
tors observe as it develops through the atmosphere. The 
photomultiplier tubes triggered by the shower define an 
arc on the sky, and, together with the position of the de- 
tector, the arc determines the so-called shower-detector 
plane. When an air shower is observed in stereo, the 
shower trajectory is in principle simply the intersection 
of the two planes. This method can be further improved 
by also taking advantage of the timing information of 
the tubes, and in our analysis the shower geometry is 
determined by a global x 2 minimization using both the 
timing and pointing information of all tubes. From mea- 
surements of laser tracks and stars in the field of view of 
the cameras we estimate that the systematic error in the 
arrival direction determination is not larger than 0.2°, 
mainly caused by uncertainties in the survey of mirror 
pointing directions. 

Various aspects of the HiRes det ector and the recon - 
struction proced ures ar e described inlBoyer et al.l l|2002j) ; 
iSadowski et"aD (|2002|k iMatthews et al.l l)2003|) . 

3. THE HIRES DATA SET 

While a ground array detector can operate year-round, 
night and day, air fluorescence detectors can only be op- 
erated on dark, moonless nights with good atmospheric 
conditions. This limits the duty cycle to about 10%. 
However, several years of observation yield a data set 
with a relatively smooth distribution in sidereal time, 
modulated by an overall seasonal variation in exposure. 

For the present analysis, we subject the HiRes stereo 
event sample to the following quality cuts. We require 
a minimum track length of 3° in each detector, an esti- 
mated angular uncertainty in both azimuth and zenith 
angle of less than 2°, and a zenith angle less than 70°. 
We additionally require an estimated energy uncertainty 
of less than 20% and x 2 /dof < 5 for both the energy and 
the geometry fit. Weather conditions which reduce the 
quality of the data are cut implicitly in the above sam- 
ple, rather than by explicit weather cuts. A total of 271 
events above 10 19 eV pass the selection criteria. A sky 
map in equatorial coordinates of the arrival directions of 
these events is shown in Figure^ 

The angular resolution of HiRes is determined using 
simulated showers. We use a full detector si mulation of 
proto n showers generated with CORSIKA 6 ijHeck et al.l 
1998) using QGSJET for the first interaction. Apply- 
ing the same cuts to the simulation data which are ap- 
plied to the real data, 68 % of all showers generated at 
10 19 eV are reconstructed within less than 0.57° of the 
true shower direction. The angular resolution depends 
weakly on energy, with the 68 % error radius growing to 
0.61° and 0.69° for showers generated at 4 • 10 19 eV and 
10 20 eV, respectively, because at higher energy, showers 
are on average farther away. The angular resolution is 
essentially constant in zenith and azimuth angle of the 
arrival direction, varying by less than 0.1°. 



Using the same simulation described above, we gener- 
ate an isotropic distribution of showers with a differential 
spectral index a = —3.0 in energy, and use the resulting 
distribution of reconstructed Monte Carlo events to de- 
termine the detector acceptance in zenith and azimuth. 
We then randomly match the local coordinates of these 
events with times during which the detector was operat- 
ing in order to generate an exposure map in equatorial 
coordinates. Figure|21 shows the distributions of the data 
and Monte Carlo events in right ascension and declina- 
tion. 

4. METHOD 

We search for small-scale clustering by performing an 
auto-correlation scan in energy and angular separation. 
Essentially, we consider the set of N events above energy 
E, count the number of pairs n p separated by less than 
9, and evaluate the probability P(N, 9) of finding this 
number or more pairs, given N and 9. We repeat this for 
a range of values for E and 9, and use the smallest prob- 
ability P m in found in the scan to identify the strongest 
clustering signal. We estimate the statistical significance 
P c h of this signal by performing identical scans over sim- 
ulated sets of isotropically distributed data, counting the 
fraction of simulated sets which yield the same or smaller 
value for P mm . 

The virtue of this approach is that by letting the en- 
ergy threshold vary, we let the scan itself determine the 
optimal balance between the better statistics of the low 
energy data set and the (presumably) smaller angular 
deflections at high energies. Furthermore, we can simul- 
taneously look for clustering both at the angular scale 
identified by AGASA and at smaller scales that take ad- 
vantage of the HiRes angular resolution. The statistical 
penalty for performing multiple searches is accounted for 
in the final evaluation of the significance P c h- 

We note that, just as in the usual two-point correla- 
tion function, higher-order multiplets are counted by the 
individual number of pairs which they contain. 

To determine the probabilities P(N, 9), we generate a 
large number of simulated data sets (typically 10 7 ) cor- 
responding to an isotropic distribution of cosmic rays. 
Specifically, we generate an event with a random ar- 
rival direction in equatorial coordinates, and accept that 
event into the simulated data set with a probability pro- 
portional to the HiRes exposure in that region of the 
sky. We then construct a table of values Pmc, where 
Pmc{N,9,u) is the fraction of data sets in which the 
first iV events contain exactly n pairs separated by less 
than 9. Then the probability P(N, 9) for observing n p 
or more pairs at (iV, 9) is simply: 



P(N,9)= ]T P M c(N,0,n) = l 



n=0 



P MC (N,e,n). 



For some combination N c and 9 C , P has a minimum: 
Pmin = P(N C ,6 C ). We identify this as the strongest 
potential clustering signal. To determine the statisti- 
cal significance, we perform the same scan over umc 
Monte Carlo data sets, finding the minimum probabil- 
ity P min = P l (Nl,9 l c ) for each trial and counting the 
number of trials n* MC for which P min < P m in- The sig- 
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Fig. 1. — Skymap (in equatorial coordinates) of the 271 HiRes stereo events above 10 19 eV examined in this study. The typical error 
radius of 0.6° is used for all events. 
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Fig. 2. — Right ascension and declination of events above 10 eV 
observed from 1999 December through 2004 January. (Data — 
points with error bars; Monte Carlo — solid line.) For right as- 
cension, x 2 /dof = 0.77; for declination, x 2 / dof = 0.73. 

nificance is finally identified as: 



- 1 ch — 7 V 7 

nivic 

that is, the chance probability of observing the value 
Pmin or less in an isotropic distribution. 

The scan is performed over the total set of N = 271 
events and over angular separations 8 from 0° to 5° in 
increments AO of 0.1°. Rather than use an arbitrary fixed 
increment AE of energy, we increase the energy threshold 
one event at a time (AN — 1). These search parameters 
were chosen a priori. While the results inevitably depend 
on t he exact choices, the dep endence is relatively small 
f see lFinlev fc W csterhoff 2004], for details and examples). 

To demonstrate the effectiveness of this method and 
the sensitivity of the HiRes detector, we apply this tech- 
nique to simulated data with clusters. First, we gen- 
erate a set of 271 events with the HiRes exposure for 
isotropic arrival directions. We then insert m pairs of 
events among the Njj highest energy events in the set to 
simulate clustering above a specific energy threshold. 

To create a pair, we pick a point in the sky for the 
source location and generate two events with arrival di- 



rections deviating from the source location according to 
a Gaussian distribution described below. These artificial 
cluster positions are chosen at random, but their dis- 
tribution is forced to reflect the overall exposure of the 
HiRes detector, so that regions with higher exposure are 
more likely to contain a cluster. The pair of events is 
then added to the original isotropic data set, replacing 
two of the original events in the set. This is repeated 
until m pairs have been inserted. The set may contain 
more than m pairs due to chance. 

For simplicity, we use a circular Gaussian distribution 
for the smearing of arrival directions around the source 
location. The width of the distribution an can be set 
equal to the angular resolution of the detector, or it can 
be set to a larger value to simulate additional smearing 
by magnetic fields. (Note that for the Gaussian distri- 
bution P(0) = (0/a 2 )e- e2 / 2a2 , the value 6 = 1.515 a 
encloses 68% of the distribution. We therefore define 
<tr = 1.515 a.) 

Table^shows the results of these simulations using the 
detector resolution (an = 0.6°), as well as three times the 
detector resolution (an = 1.8°) to simulate additional 
smearing by magnetic fields. For each choice of Nh, rn, 
and an, we generate 10 4 data sets, and scan them with 
the procedure described above to find a distribution of 
values for the significance P c h . The median and 90 th per- 
centile values of this distribution are indicated in Tabled 

The table shows, for example, that for a clustering 
signal on the an — 0.6° scale, even three pairs among 
the 47 highest energy events would typically result in 
P c h = 1.1%. The table also shows that three such pairs 
would result in P c h < 6.7% for 90% of the simulated sets. 
Thus, an actual value of P c h > 6.7% could be used to ex- 
clude the possibility that sources contributed three such 
pairs at more than the 90% confidence level. 

These results demonstrate the sensitivity to clustering 
on small angular scales. 

5. RESULTS AND DISCUSSION 

We perform the scan on the HiRes stereo sample of 
271 events above 10 19 eV. Because we start well below 
the 4 • 10 19 eV energy associated with the AGASA clus- 
tering signal, our search should safely encompass the en- 
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Table 1. Results for Simulated Clusters 



a R = 0.6° o R = 1.8' 
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a iVjy = 27, 47, and 89 events corresponds to simulated clustering 
above energy thresholds 40 EeV, 28 EeV, and 20 EeV, respectively. 
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Fig. 3. — Autocorrelation scan of the HiRes data set above 
10 19 eV. P(N, 8) is the probability of obtaining the same or greater 
number of pairs as is actually observed in the data using a maxi- 
mum separation angle 9 and searching among the N highest-energy 
events. These probabilities do not include the statistical penalty 
due to scanning. 



ergy region of interest even in the presence of a system- 
atic energy s hift of 30% between the two e xperim ents, as 
suggested bv iDe Marco. Blasi. fcOlintol 1)2003(1 . Start- 
ing at this energy does not appreciably dilute the sig- 
nificance of a clustering signal if one is found at higher 
energy, since the scan involves repeated searching with 
successively higher energy thresholds. An additional mo- 
tivation for starting at 10 19 eV is the fact that the HiRes 
angular resolution (0.6 °) is much sharper at this energy 
than AGASA's (2.8°) (iTakeda et al.lll999h . 

The results of the scan are shown in Figuredl The 
strongest clustering signal (P m in = 1-9%) is observed 
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using the energy threshold E c — 1.69 • 10 19 eV where 
we observe n p — 10 pairs separated by less than 9 C = 
2.2° within a set of N c — 120 events. The statistical 
significance of this result corresponds to P c h = 52%. 

The HiRes stereo data above 10 19 eV is therefore con- 
sistent with the null hypothesis of isotropic arrival direc- 
tions. 

Comparison with the AGASA clustering result is not 
straightforward. The HiRes stereo event sample above 
4 • 10 19 eV is still smaller than AGASA's, though how 
much smaller depends critically on the level of agree- 
ment in absolute energy scale for the two experiments. 
The possibility of a systematic energy shift of 30% 
would imply that above the rescaled energy threshold, 
(0.7) • 4 • 10 19 eV = 2.8 • 10 19 eV, HiRes has seen 47, rather 
than 27, events. More importantly, there is the question 
of how many pairs an independent data set might be ex- 
pected to contain, given the lack of an obvious source 
model and the widely varying estimates of the strength 
of the AGASA clustering. Without assuming a model 
and source strength, there is no natural way to translate 
the AGASA observation of five doublets and one triplet 
separated by less than 2.5° into a meaningful prediction 
for HiRes. 

However, what can be tested using a statistically inde- 
pendent data set is the claim that significant small-scale 
clustering is a general feature of ultrahigh energy cosmic 
ray arrival directions. The HiRes stereo data set does 
not support such a claim. We observe no statistically 
significant evidence for clustering on any angular scale 
up to 5° at any energy threshold above 10 19 eV. 

Comparing the observed value of P c h with the val- 
ues obtained from simulations in Section 4 (shown in 
Tabled, we note that if the current HiRes data above 
4 ■ 10 19 eV contained two or more pairs of events con- 
tributed by compact sources at the angular resolution 
limit of the detector, then the typical value of P c h would 
be 0.018 or less, and more than 90% of the time the value 
of P c h would be much smaller than the observed value of 
0.52. 

Results of searches for correlations with known astro- 
physical source classes will be published in a separate 
paper. 
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